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Andreev reflection in edge states of time reversal invariant Landau levels 
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We describe the conductance of a normal-superconducting junction in systems with Landau levels 
that preserve time reversal symmetry. Those Landau levels have been observed in strained hon¬ 
eycomb lattices. The current is carried along the edges in both the normal and superconducting 
regions. When the Landau levels in the normal region are half-filled, Andreev reflection is maximal 
and the conductance plateaus have a peak as a function of filling factor. The height of those peaks 
is quantized at 4 e 2 /h. The interface of the junction has Andreev edge states, which form a coher¬ 
ent superposition of electrons and holes that can carry a net valley current. We identify unique 
experimental signatures for superconductivity in time reversal invariant Landau levels. 
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Introduction. At zero temperature, Cooper pairs are 
protected against phase decoherence by time reversal 
symmetry (TRS). The most promising scenario to study 
the coexistence of superconductivity and quantum Hall 
stated requires a system that creates well separated Lan¬ 
dau levels (LLs) in the complete absence of magnetic 
flux 2 . In strained honeycomb lattices, uniform strain 
fields can couple to the electrons as a pseudo-magnetic 
field oriented in opposite directions in the two valley^E. 
As observed experimentally in graphene 5 7 and in de¬ 
formed honeycomb optical lattices^, this field can pro¬ 
duce sharp LL quantization, and at the same time, zero 
net magnetic flux at every lattice site. 

In the usual semi-classical picture, electrons moving in 
cyclotronic orbits are reflected as holes at the interface 
of a superconductor. At low energy, the holes have op¬ 
posite momentum (valley) and the same velocity of the 
incident electrons. Since particle-hole conversion changes 
the sign of the Lorentz force and effective mass, the An¬ 
dreev reflected holes and the electrons move coherently in 
the same direction of the normal-superconducting (NS) 
interface. The result is an Andreev edge state, a co¬ 
herent superposition of electrons and holes, which move 
in alternating skipping orbits along the NS interfac^^lAfll. 
Pseudo-magnetic fields nevertheless reverse their direc¬ 
tion in opposite valleys, preserving TRS. As a conse¬ 
quence, when the energy of the quasiparticles, e, is small 
compared to the Fermi energy, /it, reflected holes can re¬ 
trace the cyclotronic path of the incident electrons, and 
either form a bound state or counterpropagate along the 
same insulating edge, as shown in Fig. la. 

In this Rapid Communications, we study the trans¬ 
port across a NS junction in a honeycomb lattice with 
a discrete spectrum of TRS LLs. The current is carried 
through skipp ing orbits along the edge of the system both 
in the normal 11 and in the superconducting regions, as 
depicted in Fig. 1. We address the regime where the 
coherence length is longer than the magnetic length. In 
the limit where the energy of the quasiparticles is small 
compared to the separation of the LLs, we show that 
particle-hole conversion produces a peak in each longi¬ 
tudinal conductance plateau. Those peaks are centered 


around partial filling factors v = 4 n, n £ Z, when the 
normal LLs are half-filled, and their height is quantized 
at (n + |)8 e 2 /h. 

The Andreev edge states carry a finite charge current 
per valley along the NS interface. The valley current 
becomes asymptotically small at large energy (e /r), 
when electrons and holes move with the same group ve¬ 
locity. In the opposite regime, e -C /i, electrons and 
holes have opposite group velocities along the interface, 
and the valley current is finite. We find transport and 
spectroscopy signatures that uniquely identify proximity 
induced superconductivity in TRS LLiPI 12 ] 

Hamiltonian. In the continuum, the electronic Hamil¬ 
tonian in the presence of a pseudo magnetic-field is 13 



where v a = ^(z^o + av z ) is the projection operator into 
the valleys a = ±, and 

H a (A) = v{—iS7 + aA) • a a — n (2) 

is the Dirac Hamiltonian in each valley. a a = {a x ,aa y ) 



FIG. 1. Semi-classical picture of the edge states in TRS LLs at 
the interface with a superconducting region. Solid black lines: 
incident electrons with skipping orbits. Red lines: propagat¬ 
ing electron pairs at the edge of the superconductor; dashed 
lines: Andreev reflected holes, a) e /r regime: holes are 
retroreflected and retrace the path of the incident electrons, 
preserving their guiding center. They can form a bound state 
or counterpropagate along the same edge, b) e > fi regime: 
Andreev holes are specularly reflected and propagate along 
the NS interface as a superposition of electrons and holes. 
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is a vector of Pauli matrices, v is the Fermi velocity and 
/i is the chemical potential. 

Near the NS interface, the Bogoliubov-deGennes 
(BdG) Hamiltonian is 


'HBdG 


wher^Sl 


(H 0 (A) A(r) \ 
VA*(r) -THo^T- 1 ) 


T = 


0 a, 
a z 0 


C = v x ® o z C 


( 3 ) 

( 4 ) 


is the time reversal symmetry operator, and C the 
charge conjugation. Because the field A preserves TRS, 
T'Hq(A)T ~ 1 = 'Hq(A). In a singlet state, the elec¬ 
trons pair symmetrically across the valleys. The simplest 
Ansatz for the off-diagonal term can be written a^l 


a c>=(: o° 


A(r)v x 0 ct 0 . (5) 


Assuming a sharp NS interface, the superconductor gap 
varies abruptly at x = 0, 


A(r) 


A , for x > 0 
0 , for x < 0, 


( 6 ) 


which separates the normal (x < 0) from the supercon¬ 
ducting region (x > 0). 

In the Landau gauge, A = (—By, 0), the Hamiltonian 
in the valleys can be written as 

V 

Ua = v^£~ a ( idi ° v + ^ ( 7 ) 

where £ = tBk — y / £b is a dimensionless coordinate with 
guiding center X a = isk, and £b — \fh/Be is the mag¬ 
netic length (restoring H). By convention, we define the 
momentum along the edge for each valley as k x = ak, 
with k > 0. The electronic wavefunction in the normal 
region moving towards (away from) the interface in val¬ 
ley a = + (—) takes the form r ip a (x,y) = e mkx 'I> a (y), 
where 'F a (?/) = ( (j>A,a , <pBa) is a two component spinor in 
sublattice A and B of the honeycomb lattice. 

Edge states. In order describe the edge states, we 
introduce a mass term potential M(y)v o (g) a z , where 
M(y ) = W for y < 0 and M = 0 for y > 0. In the limit 
W —> oo, this potential describes the edge of the system 
at y = 0. The electrons will move in skipping orbits along 
the y = 0 line under the influence of a pseudo-magnetic 
field. In the normal side of the NS interface, 


[H a + M(y)a z \^ a = e'i’a, (8) 


where 'L ct (y) is a two component spinor in the sublattice 
basis. Multiplying Eq. Q on the left by the charge 
conjugated form of T~L a + M(y)a z , where y —> —y, this 
equation assumes a diagonal form, 


-A-7 ) e-M\y) + (s + yf--c 


^a = 0. (9) 



FIG. 2. Energy spectrum versus guiding center k x lB at the 
edge (y = 0) for y = 1.5 (v = 4) and A = 0.5. Energy 
scales in units of \/2 v/Ib- a) Normal region (x < 0). b) 
superconducting region (a; > 0). 


For y > 0, where M = 0, the energy spectrum of the LLs 
at the edge is 

e(k x ) = s(V2v/£ B )V\ n {k x )\ - y (10) 


where n(k x ) is a real number and s = sign[n(£; x )]. The 
corresponding eigenvectors in the two valleys are of the 
form 


*a,n(0 


D n (k x )—i 



( 11 ) 


where D n ^ k ^(x) are parabolic cylinder functions. The 
determination of n(k x ) follows from enforcing boundary 
conditions at the edge and results in a discrete number 
of edge states shown in Fig. 2a. For definiteness, we con¬ 
sider the case of a zigzag edge, where 0 b,+ (^_b/V 2) = 
0, although similar conclusions apply to any choice of 
boundary conditions. 

In the superconducting region, we can decompose the 
BdG Hamiltonian © into two identical blocks of 4 x 4 
matrices, where 


( H a 

\A*(7 0 


Act 0 \ 

-Ha) 


<F a = .E-E 


( 12 ) 


is the reduced BdG equation and = ('ft e ,a'&h,a) is a 
4 component spinor including electron and hole states. 
At y > 0, where the wavefunctions are finite (M = 0), 
the solution follows from squaring (12) with the charge 


conjugated BdG Hamiltonian, which results in the differ¬ 
ential equation 


1 d 2 1 f 2 

2 « z 


2 yM 


= 0 , 


(13) 


where 


M = 


e - If t z -Atj 0 

-Afj 0 —e - \cr z 


is a 4x4 matrix in the particle-hole basis, and A is as¬ 
sumed real. The four component spinor that satisfies Eq. 


(13) and hence (12) is 
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<f> 


ot. 


Pc& a ,n.(Z) \ 

/3_ a 'F a ,„ fl (£) J 


(14) 


where f3 a = \(\ + ct\Jl — A 2 /E 2 ), with a = ± index¬ 

ing the valleys. The energy spectrum in the supercon¬ 
ducting edge is given by 


E(k x ) = \j(^s{V2v/£ B )A\n s {kx)\~ +A 2 , (15) 

where s = sign[n s (fc x )] and n s (k x ) is a real number to 
be found from the boundary condition at the edge (y = 
0), in superconducting side. Imposing similar boundary 
conditions <j) B ,+ {k£ B /V 2) = 0, the energy spectrum at 
the superconducting edge is shown in Fig. 2b. 

Transport across the NS junction. In the normal re¬ 
gion, the electron and hole like excitations are decou¬ 
pled. Near the NS interface, the normal edge state can 
be written as a superposition of the wavefunctions of the 
incident electron, ipj(x,y), the reflected electron in the 
opposite valley, ipg(x, y), and the Andreev reflected hole, 
tp//{x,y). For simplicity, we assume from now on all 
length scales to be in units of the magnetic length £ B 
and all energy scales to be in units of y / 2v/l B . 

The largest contribution to scattering comes from 
states at integer values of n(k) and n s (k), where the 
density of states is the largest. In the four component 
particle-hole basis, the electron wavefunctions are 


Ki^y)=^ + ’ nA O i (16) 


and 


Q 

v) = E r h ( o j ' y] ) e ~ ikjX ’ ( 17 ) 


where the nearest integer function n e = nint[(e + p,) 2 ] £ 
Z sets the index of the highest occupied band, with C Ue = 
2 n e + 1 the number of channels crossing the Fermi level 
at the edge. The wavefunction of the Andreev reflected 
hole is 


^hA x ’V) = 


Cn h 

E r tj 


^-,n h (kj - y) 


—iki 


(18) 


where nh = nint[(e — y) 2 ] £ Z, with C nh equivalently 
the number of edge channels for hole states. 

In the superconducting region, the electron-like and 
hole-like solutions can be written as, 


r s ,i 


c„, 


E A h 


3 


/3a’F an+ (kj y) \ nikjX 

/3- a ^ a ,n+ {kj -y) ) 


(19) 




FIG. 3. Top panel: Andreev reflected hole probability |r,t| 2 
per channel as a function of the normal filling factor v. The 
superconductor gap A (in units of \/2v/£b) ranges from 0.01 
(red line) to 0.2 in steps of 0.025. Bottom: corresponding 
longitudinal conductance G at the edge as a function of v. 
The solid line plateaus: normal conductance (A = 0). Dotted 
lines: upper bound for the conductance, which is quantized 
at 4 e 2 /h. In all curves, the temperature T = 0.02(\/2 v/£b)- 


where n a = mnt[{\/E 2 — A 2 + ay) 2 ] £ Z. The reflected 
electron and hole amplitudes can be calculated by match¬ 
ing the amplitudes at the interface (x = 0), 

^t+^e +V>h =i’s+^S- ( 20 ) 


In the limit of large coherence length, £ = v/A £ B , 
and e <C v/i B , we can neglect scattering processes be¬ 
tween different modes. Also, the number of edge chan¬ 
nels is the same for electrons and holes in the two sides 
of the junction, n K = nh = n a . In this regime, the 
solution of Eq. (201 is rf = 0 and rf = /?+//?_ = 


E/A — A A 2 /A 2 — 1. When E < A, rf is complex and 
the total amplitude of the Andreev reflection is |r) 4 | 2 = 1. 

In fig. 3, we show the amplitude of the Andreev re¬ 
flection versus the filling factor of the normal LLs. In 
the normal region, when the LLs are well separated, the 
Fermi distribution of the quasiparticles is equal to the fill¬ 
ing fraction of the highest occupied LL, (e e / T + l) -1 = 
f{v) = (i//4 + |)mod(l) £ [0,1], where T is the temper¬ 
ature. The Andreev amplitude is rf = 0(z/), where 


0M = 


ME 

A 



( 21 ) 


with 


eM=T AAw))’ (22) 

From the Blonder-Tinkham-Klapwijk formula 1 ^, the 
differential conductance at the NS junction can be writ- 
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FIG. 4. Energy spectrum e along the NS interface (x = 0) 
versus guiding center kyls . Negative guiding centers ( k y (.B < 
0) describe states in the normal region. k y £s > 0 correspond 
to states in the superconducting one. Left curves: normal 
edge states. Center right curves: low energy Andreev edge 
states, a) /x = 0 [u = 0), e > g regime for finite A. Electron 
and hole states have the same group velocity v g = de/d k y and 
propagate in the same direction of the interface; b) fj/A = 2: 
electrons and Andreev reflected holes have the same guiding 
center and opposite group velocities along the interface for 
e -C /x, forming an Andreev bound state. 


ten as 


G = x Dh^I'+i^i 2 ) * x ( n + IY 1 + |0(ty)|2] ' 

i —1 ' ' 

( 23 ) 

The peaks in the differential conductance are shown in 
the bottom panel of Fig. 3 as a function of the filling 
factor of the normal region v. The solid line plateaus 
represent the conductance in the absence of Andreev re¬ 
flections. The different curves show the conductance for 
different values of the normalized gap A ranging from 
0.01 (red line) to 0.2 in 0.025 steps. The peaks appear 
at partial fillings v = An, n £ Z, and their height is 
quantized at ( 2 n+ l)4e 2 /h. At those fillings, the normal 
LLs are particle-hole symmetric and Andreev reflection 
is maximal. At integer fillings v = 4 (n + |), when the 
Fermi level is in the middle of the LL gap, Andreev re¬ 
flection is suppressed and the conductance is quantized 
by half, at (2 n + 1)2 e 2 /h. 

Andreev edge states. More insight on Andreev reflec¬ 
tion and the electronic states near the interface can be 
obtained by considering the current flowing parallel to 
the interface. In the configuration where the insulating 
edge is zigzag, the NS interface has armchair character. 
Its states are formed by a superposition of eigenstates on 
both valleys. Using the Landau gauge A = (0 ,Bx), the 
wavefunction in the normal side of the interface can be 
written as 


V’H (x,y) = 


^e,a^Q,n e (e) ( x 
a,n,h(e) ( x 


ky) 

ky) 


e aik v y 


(24) 


where n e (e) = (e + /i ) 2 and n/,(e) = (e — /x ) 2 are real 



FIG. 5. Andreev edge states at the NS interface. Solid lines: 
electron cyclotronic orbits; dashed: Andreev reflected holes. 

a) e/fi —> 0 regime: electrons and holes form a bound state. 

b) Intermediate regime, e < min(^t, A): electrons are retrore- 
flected into holes with group velocity v g = de/dk y having 
opposite sign, c) e > /x regime: electrons are specularly re¬ 
flected into holes, which propagate along the same direction. 


numbers. In the superconducting side, 


Tps\\(x,y) = Y, A i'y 


/^7^a,n a (e) ( x ky) \ ^ aik y y 

ft—a,n a (e) ( x ky) J 


(25) 

with 7 = +(—) indexing electron(hole)-like states, and 

n a( £ ) = (Ve 2 — A 2 + an) . Matching the amplitudes 
V>||( 0 , y) = i/jsw ( 0 , y) and the derivatives d x ip\\(x,y) = 
d x ips\\ ( x i y) at x = 0 yields eight linear equations for 
eight unknown coefficients V % . with i = 1,..., 8 . This set 
of equations can be expressed in matrix form as Q- V = 0. 
Non-trivial solutions require that Det(Q) = 0. From this 
condition, we numerically extract the spectrum of exci¬ 
tations e(k y ) near the NS interface, as shown in Fig. 4. 

The interface states are an admixture of two type of 
modes: i) normal edge states formed by conventional 
skipping orbits moving in one direction, and ii) Andreev 
states, which are coherent superpositions of particles and 
holes. The first mode is connected to bulk LL energies in 
the normal region as k y —> — oo. The second one has no 
correspondence with the bulk LLs in the normal region, 
and appears around kyls ~ 0 and also inside the super¬ 
conducting region, for positive guiding centers (k y > 0 ). 

In panel 4a, we plot the interface modes at /x = 0 
(u = 0), for finite A. In this regime, e > y = 0, the group 
velocity of the Andreev edge states v g = de/dk y > 0 for 
both electrons and holes (center ri ght curves), which are 
specularly reflected at the interface 16 - and move in alter¬ 
nating skipping orbits in the same direction (see fig. 5c). 
In the limit e /x , electrons and holes have the same 
velocity and guiding center, and carry zero net charge per 
valley. In all other cases, their velocities are different, re¬ 
sulting in a net valley current along the NS interface. In 
panel 4b, we plot the energy of the modes for fi/A = 2. 
In the regime e < min(/x, A) we find numerically that the 
Andreev reflected holes and electrons have group veloci¬ 
ties v g with opposite signs (fig. 5b). In the limit e//x —> 0, 
the holes retrace the path of the incident electrons, form¬ 
ing an Andreev bound state schematically shown in fig. 


5a. 


In summary, we have derived transport and spec¬ 
troscopy signatures of proximity induced superconductiv¬ 
ity in TRS LLs. We found the longitudinal conductance 
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as a function of the filling factor across an NS junction, 
and showed that it is quantized at (2 n + l)4e 2 //i in the 
n-th LL at half-filling, when Andreev reflection is maxi¬ 
mal. We also showed that the NS interface has Andreev 
edge states, with unique spectroscopic features. 
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